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Abstract. It is well-known that solutions to the conformal formulation of the Ein- 
stein constraint equations are unique in the cases of constant mean curvature (CMC) and 
near constant mean curvature (near-CMC). However, the new far-from-constant mean 
curvature (far-from-CMC) existence results due to Hoist, Nagy, and Tsogtgerel in 2008, 
to Maxwell in 2009, and to Dahl, Gicquaud and Humbert in 2010, are based on degree 
' theory rather than on the (uniqueness-providing) contraction arguments that had been 

fN| . used for all non-CMC existence results prior to 2008. In fact, Maxwell demonstrated 

q " in 2011 that solutions are non-unique in the far-from-CMC case for certain types of 

D . low-regularity mean curvature. In this article, we investigate uniqueness properties of 

' solutions to the Einstein constraint equations on closed manifolds using tools from bi- 

furcation theory. For positive, constant scalar curvature and constant mean curvature, 
we first demonstrate existence of a critical energy density for the Hamiltonian constraint 
with unsealed matter sources. We then show that for this choice of energy density, the 
linearization of the elliptic system develops a one-dimensional kernel in both the CMC 
and non-CMC (near and far) cases. Using Liapunov-Schmidt reduction and standard 
tools from nonlinear analysis, we demonstrate that solutions to the conformal formula- 
tion with unsealed data are non-unique by determining an explicit solution curve, and by 
analyzing its behavior in the neighborhood of a particular solution. 
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1. Introduction 

In this paper we demonstrate that solutions to the Einstein constraint equations on a 
3 -dimensional closed manifold (M.,g a b) with no conformal killing field are non-unique. 
More specifically, we show that solutions to the conformal formulation of the constraint 
equations with an unsealed matter source on (A4, g) exhibit non-uniqueness in the case 
that the scalar curvature is positive and constant. Letting k ab be a (0, 2) tensor and R and 
D be the scalar curvature and connection associated with g ab , the constraint equations 
take the form 

R + k 2 - k ab k ab = 2kp, (1.1) 
D a k + D b k ab + Kf = 0. (1.2) 

Equation (11.11) is known as the Hamiltonian Constraint and (11.21) is known as the mo- 
mentum constraint. 

Equations (|l.ll) and (11.21) form a system of coupled elliptic partial differential equa- 
tions. When one attempts to solve the constraint equations they are faced with the prob- 
lem of having twelve pieces of initial data and only four constraints. One solution to 
this problem is to attempt to parametrize solutions to (|l.ll) and (11.21) by formulating the 
constraints so that eight pieces of initial data are freely specifiable while four are deter- 
mined by dl.lh - dl.2h . The conformal transverse traceless (CTT) decomposition and the 
conformal thin sandwich method (CTS method) are standard ways of doing this. The 
extended conformal thin sandwich method (XCTS method ) is popular among numerical 
relativists and reformulates (11.11) and (II .21) as a coupled system of 5 elliptic equations. 
In the CTT method one decomposes k ab into its trace or mean curvature and trace free 
part and then scales this trace free tensor, the metric g ab and the source terms p and j by 
judicious choices of some power of a positive, smooth function </>. The choice of scaling 
power for each term is typically made to simplify the analysis of the resulting system. 
In particular, one chooses powers to eliminate terms involving (D a (f))/(f) and so that the 
system decouples when the mean curvature is constant. 

It is well-known that solutions to the CTT formulation of the constraint equations with 
scaled data sources are unique in the event that the mean curvature is constant (known as 
the "CMC case"), or near constant (the "near-CMC case"); cf. 011 [II [BE CGI. Prior 
to 2008, all non-CMC existence results were only possibly for the near-CMC case, and 
were established using contraction arguments, which provided uniqueness for free once 
existence was established. However, beginning in 2008 with the first true "far-from- 
CMC" (the non-CMC case without near-CMC restrictions) existence result in ||9), all 
far-from-CMC results to date flUQISEEH are based on a variation of the Schauder Fixed- 
Point Theorem. This also includes the more recent work Q, which uses the Schauder 
framework from JHELOlEEH as part of a pseudo-variational argument. As a result, little is 
known about uniqueness of far-from-CMC solutions. In fact, in 2011 Maxwell demon- 
strated in |[T5ll that solutions of the CTT formulation of the constraint equations are non- 
unique in the far-from-CMC case for certain families of low regularity mean curvatures. 
However, as noted by Maxwell in 03), the discontinuous mean curvature functions con- 
sidered by Maxwell in 031 fall outside of the best existing non-CMC rough solution 
theory established in IfTOl . 

In 0/7), Pfeiffer and York provided numerical evidence for non-uniqueness of the 
XCTS method on an asymptotically Euclidean manifold. In 0), Baumgarte, O'Murchadha, 
and Pfeiffer conjectured that the non-uniqueness demonstrated by Pfeiffer and York was 
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related to the fact that certain terms in the momentum constraint related to the lapse func- 
tion have the "wrong sign", which prevents an application of the maximum principle. To 
support their claim, the authors of JU analyzed a simplified system corresponding to 
a spherically symmetric constant density star and explicitly constructed two branches 
of solutions. In their analysis they proved that solutions to the Hamiltonian constraint 
(11.11) with an unsealed matter source are non-unique. Then in [fT9l , Walsh generalized 
the work in [4j by applying a Liapunov-Schmidt reduction to both the Hamiltonian con- 
straint with an unsealed matter source and to the XCTS system on an asymptotically 
Euclidean manifold. However, Walsh relied on the assumption of the existence of a criti- 
cal density for which the linearization of these two systems developed a one-dimensional 
kernel. Here we extend the work of Walsh by applying a Liapunov-Schmidt reduction 
to the CTT formulation of the constraint equations on a closed manifold. We explicitly 
construct a critical, constant density in the event that the scalar curvature is positive and 
constant and the transverse traceless tensor has constant magnitude. For this particular 
density, we then show that solutions to the CTT formulation with an unsealed density are 
non-unique. 

As in BH [191, we consider a less standard conformal formulation of the constraints 
by allowing unsealed matter sources p and j. However, as opposed to considering the 
CTS and XCTS formulations as in ffTTl H Q3D, we consider the CTT formulation. By 
decomposing our initial data 

Kb = Lb + 7^g a bt, (1.3) 

where f = k a b9 ah is the trace and l a \> is the traceless part, making the following conformal 
rescaling 

9ab = <p 4 9ab, L = <p- W l ab , t = T, (1.4) 

and then decomposing 

Lb = {°ab + (£w) ab ), (1.5) 

where D a a ab = and 

(Cw) ab = D a w b + D b w a - -{D c w c )g ab 

3 

is the conformal Killing operator, we obtain the following unsealed conformal refor- 
mulation of (|l.ll) and (11.21) that we will analyze 

-A<P+Ir<P + -3-r 2 5 - \{a ab + (Cw) ab )(a ab + (£w) a6 )0- 7 - 2tt P 5 = 0, (1.6) 

o iZ o 

- D b (Cw) ab + -D a r<f + /tj a 10 = 0. 
3 

Our non-uniqueness results for (|1.6I) are of interest for a number or reasons. Most 
immediately, our analysis shows that the formulation (11.61) is unfavorable due to the non- 
uniqueness of solutions. Therefore, for a given system, if the CTT formulation with a 
scaled matter source leads to a set of constraints that is suitable for analysis, which it 
usually does, then one should use the scaled formulation. However, it is not always the 
case that the conformal formulation with scaled matter sources is the ideal formulation 
for a given source. In the case of the Einstein- scalar field system, the conformal formu- 
lation that is most amenable to analysis takes on a form very similar to the system (11.61) 
[[5]|. In addition, it is the hope of the authors that these results will provide additional 
insight into the non-uniqueness phenomena associated with the CTT formulation in the 
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far-from-CMC case [fl"5ll and with the non-uniqueness phenomena analyzed by Pfeiffer 
and York [171, by Walsh [ 19] and by Baumgarte, O'Murchadha, and Pfeiffer H. In par- 
ticular, the analysis conducted in this article clearly demonstrates the effect that terms 
with the "wrong sign", as discussed in [4], have on the non-uniqueness of the conformal 
formulations of the constraints. In the case of (|1.6I) . the negative sign in front of the term 
27rp0 5 is undesirable given that it prevents the semilinear portion of the Hamiltonian 
constraint from being monotone and the corresponding energy from being convex. By a 
maximum principle argument, we will see in section 147X1 that it is this term that directly 
contributes to the non-uniqueness properties of ( 11.6b . 

The rest of this paper is organized as follows. In section [2] we introduce the func- 
tion spaces that we will use and some basic concepts from functional analysis. Then 
we discuss the Liapunov-Schmidt reduction that we use to prove non-uniqueness. The 
statements of the main results of this paper can be found in section [73] The remainder 
of this paper is then devoted to proving these results. The foundation for our argument 
is developed in sections 14.11 and 14.21 In section 14.11 we demonstrate the existence of a 
critical, constant density p c such that if g a f, has positive, constant scalar curvature, |er| 
is constant and j" = 0, the Hamiltonian constraint in (11.6b will have a positive solution 
if P < Pc an d will have no positive solution if p > p c . Then in section [4721 we use the 
properties of p c to show that there exists a function <\> c at which the linearizations of the 
uncoupled Hamiltonian operator (CMC case) and coupled system (non-CMC case) have 
one-dimensional kernels. The existence of a one-dimensional kernel then allows us to 
apply the Liapunov-Schmidt reduction in section [5711 in the CMC case and in section 15721 
in the non-CMC case. In particular, in section 15.11 we determine an explicit solution 
curve for (11.6b that goes through the point (0 C , 0) in the CMC case. An analysis of this 
curve then implies the non-uniqueness of solutions to (|1.6b when the mean curvature is 
constant. Similarly, in section 15.21 we also determine an explicit solution curve for the 
full, uncoupled system (11.6b through a point of the form ((<f> c , 0), 0). Again, an analysis 
of this curve reveals non-uniqueness in the event that the mean curvature is non-constant. 

2. Preliminary Material 

In this section we give a brief definition of the function spaces, norms and notation that 
we will use in this article and then discuss some basic concepts from functional analysis 
and bifurcation theory that will be necessary going forward. 

2.1. Banach Spaces, Hilbert Spaces and Direct Sums. We introduce the fundamental 
properties of the function spaces with which we will be working. We will primarily be 
working with Banach spaces, however at times we will need to consider these spaces as 
subspaces of a Hilbert space. For convenience, we present the basic definitions of these 
general spaces and define the direct sum of two vector spaces, which will be necessary 
in our non-uniqueness analysis. 

The basic space that we will be working with is a Banach space, where a Banach 
space X is a complete, normed vector space. If the norm || ■ || on X is induced by an 
inner product, we say that X is a Hilbert Space. One can form new Banach spaces and 
Hilbert spaces from preexisting spaces by considering the direct sum. 

Definition 2.1. Suppose that Xi andX 2 are Banach spaces with norms \\-\\x 1 and || • ||x 2 - 
Then the direct sum X\ © X 2 is the vector space of ordered pairs (x, y) where x G X\, 
y G Xi and addition and scalar multiplication are carried out componentwise. 

We have the following proposition: 
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Proposition 2.2. The vector space X\ © X 2 is a Banach space when given the norm 

\\(x,y)\\ Xl(BX2 = (\\x\\ 2 Xl + \\y\\ 2 X2 )K (2.1) 



Proof. This follows from the fact that || ■ and || ■ \\ x , 2 are norms and the spaces X\ 
and X 2 are complete with respect to these norms. □ 

We have a similar proposition for Hilbert spaces. 

Proposition 2.3. Suppose that Hi andH 2 are Hilbert spaces with inner products (-, 
and (•, -)u2- Then the direct sum Hi © H 2 is a Hilbert space with inner product 

((w, x), (y, z)) Hl(Bn2 = (w, y) Hl + (x, z) n . 2 . (2.2) 

Proof. That (-, -} nien2 is an inner product follows from the fact that (■, and (•, -) H2 
are inner products. The expression 

\\(u,v),(u,v)\\n l ®u a = v 7 ((u,v), (u,v)) ni ®H2, 
is a norm on Hi © H 2 that coincides with the norm in Proposition 12.21 in the event that 
the norms on X\ and X 2 are induced by inner products. □ 

See [20] for a more complete discussion about the direct sums of Banach spaces. 



2.2. Function Spaces. Let E denote a given vector bundle over M.. In this paper we 
will consider the Sobolev spaces W k ' p (E), the space of /c-differentiable sections C k (E), 
and the Holder spaces C k ' a (E) where k G N, p > 1, a G (0, 1) and E will either be 
the vector bundle M. x R of scalar- valued functions or T s r Ai, the space of (r, s) tensors. 
Note that all of these spaces with the following norm definitions are Banach spaces and 
the space W k ' 2 (E) is a Hilbert space for fceN. 

Fix a smooth background metric g a b and let f ^ ... be a tensor of type r + s. Then at 
a given point x E M, we define its magnitude to be 

\v\ = (v ai '-' bs v ai> ... jbs )^, (2.3) 

where the indices of v are raised and lowered with respect to g a f,. We then define the 
Banach space of fc-differentiable functions C k (Ai x R) with norm || • ||^ to be those 
functions u satisfying 

k 

= sup \IPu\ < oo, 
i=0 xeM 

where D is the covariant derivative associated with g a b. Similarly, we define the space 
C k (TgA4) of fc-times differentiable (r, s) tensor fields to be those tensors v satisfying 

\\v\\k < oo. 

Given two points x,y G Ai, we define d(x,y) to be the geodesic distance between 
them. Let a E (0, 1). Then we may define the C°' a Holder seminorm for a scalar- valued 
function u to be 

Using parallel transport, this definition can be extended to (r, s)-tensors t> to obtain the 
C fc,a seminorm [u)k, a @. This leads us to the following definition of the C fc ' Q (A / l x R) 
Holder norm 

\\u\\k,a = \\u\\ k + [u] k , a 

for scalar- valued functions, and we may define the C k ' a (T s r Ai) Holder norm for (r, s) 
tensors in a similar fashion. 
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Finally, we will also make use of the Sobolev spaces W k ' p (Ai x R) and 
W k ' p (TgA4) where we assume k e N and p > 1. If dV g denotes the volume form 
associated with g a b, then the L p norm of an (r, s) tensor is defined to be 



M 



v\ p dV g \ . (2.4) 



We can then define the Banach space W k ' p (M x R) (resp. W k ' p (Tj M)) to be those 
functions (resp. (r, s) tensors) v satisfying 



k 



1'%*= i X) lip 1 " x - 



0=0 

The above norms are independent of the background metric chosen. Indeed, given 
any two metrics g a b and g a b, one can show that the norms induced by the two metrics 
are equivalent. For example, if D and D are the derivatives induced by g a b and g ab 
respectively, then there exist constants C\ and C 2 such that 

Ci|MU,s — ll M IU,9 — ^alpllfe,^, 

where || • ||fe :S denotes the C k (Ai) norm with respect to g. This holds for the W k ' p and 
C k ' a norms as well. We also note that the above norms are related through the Sobolev 
embedding theorem. In particular, the spaces C k,a and W l,p are related in the sense that 
if n is the dimension of M. and u E W l ' p and 

n 

k + a < I , 

V 

then u G C k ' a . See [2] El [81 [161 for a complete discussion of the Sobolev embedding 
Theorem, Banach spaces on manifolds, and the above norms. 

2.3. Adjoints and Projection Operators. Solutions to the coupled system (11.61) satisfy 

F(x,w) = 0, (2.5) 

where F:XxF->Zisa nonlinear operator between Banach spaces. This allows us 
to use basic tools from functional analysis to analyze our problem. In particular, we will 
repeatedly need to consider the linearization D x F(x, w), its adjoint, and projections onto 
subspaces determined by these operators. Later on in the section when we introduce the 
Liapunov-Schmidt reduction, we will use the kernel of the linearization D x F(x , w ) at a 
point (xo, Wo), the kernel of the adjoint, and projection operators onto these subspaces, to 
decompose X and Y in a manner that will greatly simplify our analysis. Here we briefly 
discuss the adjoint and projection operators. See EOl for a more complete discussion of 
these topics and see Appendix [7J] for a discussion of Frechet derivatives. 

2.3. 1 . The Adjoint and Properties. Suppose that H is a Hilbert space with inner product 
(-,•}. Then if A : H — > H is a linear operator, the Riesz Representation Theorem implies 
that there exists a unique operator A* that satisfies 

(Ax, y) = (x, A*y) for all x,yeH. (2.6) 

If R(A) denotes the range of A and ker(A) denotes the kernel, then the operator A* 
satisfies the following properties: 

1) ker(A*) = R(A) L (2.7) 

2) (ker(A*)) L = R[A). (2.8) 
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2.3.2. Projection Operators and Fredholm Operators. Now assume that X C H is a 
Banach space contained in a Hilbert space Given a subspace t^cl, the projection 
P onto is a bounded linear operator P : X — > V that satisfies P 2 = P. In particular, 
if V is a finite-dimensional subspace spanned by the orthonormal basis v\, - ■ ■ ,v n , then 
we can easily construct the projection onto V by the formula 

n 

Pu:J2( u ,Vi)vi, (2.9) 

i=l 

wherew G Xand(-,-)is the inner product on H. Note that Pis just the normal projection 
operator from % to V restricted to X. 

We end the section by introducing one more definition that will be important in the 
following section. A Fredholm operator is a bounded linear operator A : X — > Y 
where X and Y are Banach spaces such that dimker(A) and dimker(A*) are finite- 
dimensional and R(A) is closed. Given a nonlinear operator F : U — > Y where U C X, 
we say that F is a nonlinear Fredholm operator if it is Frechet differentiable on U and 
D x F{x) is a Fredholm operator. 

Notice that if A is a Fredholm operator, then ker(yl*)- L = R(A) and furthermore, 
the fact that ker(A) and ker(A*) are finite dimensional allows one to define projection 
operators P and Q onto ker(A) and ker(A*) to decompose X and Y. As we will see, 
these properties make Fredholm operators ideal candidates for bifurcation analysis. 

2.4. Elements of Bifurcation Theory. We now present some basic concepts from bi- 
furcation theory that will be essential in obtaining our non-uniqueness results. In partic- 
ular, we give a formal definition of a bifurcation point and then present the Liapunov- 
Schmidt reduction. This reduction allows one to reduce a nonlinear problem between 
infinite-dimensional Banach spaces to a finite-dimensional or even scalar-valued prob- 
lem. Therefore it greatly simplifies the analysis and will serve as a basic tool for us 
going forward. The following treatment is taken from [fT3l and H. 

Suppose that F : U x V — > Z is a mapping with open sets U C X, V C A, where X 
and Z are Banach spaces and A = K. We let x E X and A G A. Additionally assume 
that F(x, A) is Frechet differentiable with respect to x and A on U x V. We are interested 
in solutions to the nonlinear problem 

F(x,X) = 0. (2.10) 

A solution of (12.101 ) is a point (x, A) G X x A such that (12.101 ) is satisfied. 

Definition 2.4. Suppose that (x ? Ao) is a solution to (12.101) . We say that \ is a bifur- 
cation point if for any neighborhood U of (xo, Ao) there exists a A G A and X\,x% G X, 
x\ ^ X2 such that (xi, A), (x2, A) G U and (xi, A) and (x-2, A) are both solutions to 

dm . 

Given a solution (xo, Ao) to (12.101) . we are interested in analyzing solutions to (12.101) 
in a neighborhood of (xo, Ao) to determine whether it is a bifurcation point. One of the 
most useful tools for this is the Implicit Function Theorem 17 .51 This theorem asserts that 
if D x F(x , A ) is invertible, then there exists a neighborhood U\ x V\ C U x V and a 
continuous function / : V\ — > U\ such that all solutions to (12.101) in U\ x V\ are of the 
form (/(A), A). Therefore in order for a bifurcation to occur at (x , A), it follows that 
D x F(xq, Aq) must not be invertible. 
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2.4.1. Liapunov-Schmidt Reduction. The following discussion is taken from [fT3l . Let 
X, A and Z be Banach spaces and assume that U C X, V C A. For A = A , we require 
that the mapping F : U x V — > Z be a nonlinear Fredholm operator with respect to x; 
i.e. the linearization D X F(-, Ao) of F(-, Ao) : U — > Z is a Fredholm operator. Assume 
that F also satisfies the following assumptions: 

F(x ,A )=0 for some (x , A ) G U x V, (2.11) 
dimker(D x F(x , A )) = dimker(D x .F(x , A )*) = 1. 

Given that D x F(xq, Ao) has a one-dimensional kernel, there exists a projection oper- 
ator P : X — > X\ = ker(D x F(x , Ao)). Similarly, one has the projection operator 
Q : Y — )• Y 2 = ker(D x F(x , A )*). This allows us to decompose X = X x © X 2 and 
Y = Y X @Y 2 where Y"i = R(D x F(x , A )). We will refer to the decomposition X 1 © X 2 
and Yi © F 2 induced by D x F(x , A ) as the Liapunov decomposition, and we see that 
F(x, A) = if and only if the following two equations are satisfied 

QF(x,\) = 0, (2.12) 
(I -Q)F(x,\) = 0. 

For any x E X, we can write x = v + w, where v = Px and w — (I — P)x. Define 
G : U x x W 1 x V x ->■ Y x by 

u>, A) = (/ - Q)F(v + w,X), where (2.13) 
Z7i C Xi, C X 2 , ViCR and 
v = Px eU u w = (I - P)x e W u 

and Ux, W\ are neighborhoods such that Ui + W\<zU C X. 

Then the definition of G(v, w, A) implies that G(v , w , A ) = and our choice of 
function spaces ensures that 

D w G(v , wo, A ) — {I — Q)D x F(x , A ) : X 2 -> F l5 

is bijective. The Implicit Function Theorem then implies that there exist neighborhoods 
U 2 C U\, W 2 C Wj. and C Vj. and a continuous function 

ip : U 2 x V 2 ^ W 2 such that all solutions to G(v, w, A) = (2.14) 
in J7 2 x VF 2 x V2 are of the form G(v, ip(v, A), A) = 0. 

Insertion of the function ip(v,X) into the second equation of (12.121) yields a finite-dimensional 
problem 

$(v,\) =QF(v + 4>(v,\),\) = 0. (2.15) 

We observe that finding solutions (d, A) to (12.151) is equivalent to finding solutions to 
F(x,\) = in a neighborhood of (x ,X ). We will refer to the finite-dimensional 
problem (12.151 ) as the Liapunov-Schmidt reduction of (12.10I ). 

Given that ker(D x F(x , A )) is spanned by v , then we can write v = sv + v . Sub- 
stituting this into (12.151 ) we obtain 

$(s, A) = QF{sv + v + tfj{sv + v , A), A) = 0. (2.16) 

Using the reduction (12.161) and another application of the Implicit Function Theorem, one 
obtains the following theorem taken from |fT3l , which allows us to determine a unique 
solution curve through the point (xq, Ao). We also include the proof for completeness. 



NON-UNIQUENESS OF SOLUTIONS TO THE CONFORMAL FORMULATION 9 

Theorem 2.5. Assume F : U xV — >■ Z is continuously differentiable onU x V C X x R 
and that assumptions (12.1 II ) hold. Additionally we assume that 

D x F(x , A ) i R(D x F(x , A )). (2.17) 

Then there is a continuously differentiable curve through (x , A ); that is, there exists 

{(x(s), X(s)) | s e (-6, 6), (x(0), A(0)) = (a; , A )}, (2.18) 

such that 

F(x{s),\(s)) = for se (-6,6), (2.19) 

and all solutions ofF(x, A) = in a neighborhood of (x , A ) belong to the curve (12.181 ). 

Proof. Let x = v + w = v + ip(v , A ). Differentiating (|2.15l) with respect to A we 
obtain 

D A <&KAo)= (2.20) 
QD x F{x Q , \o)D x i>(v , A ) + QD x F{x , A ) = QD x F{x , A ) ^ 0, 

where (12.201) is nonzero due to the extra assumption (12.171) . The above expression sim- 
plifies due to the fact that that 

D x F(x ,X )D x i)(v ,X ) e R(D x F(x ,X )), 

and Q is the projection onto ker(D x F(x , A )*). 

The fact that D x &(v , A ) ^ and that X 1; Y 2 and R are one-dimensional implies that 
we may apply the Implicit Function Theorem to $(t>, A) to conclude that there exists a 
continuously differentiable 7 : U 2 — > V 2 C K. such that 

7(v ) = A and $O,7(¥))=0 for all veU 2 dX l . (2.21) 
Therefore our reduced equation (12.151) becomes 

7(1;)) = QF(t; + 7(«)), 7(«)) = 0, (2.22) 
where solutions to (|2.22l) are of the form 

x(v) = v + ip(v , 7(f)) and A(i>) = 7(f). (2.23) 

By writing t> = sv + vo as in (|2.16l) and inserting this into (12.231) . we obtain our solution 
curve 

x(s) =v + sv + ip(v + sv , 7(^0 + sv )), (2.24) 
\(s) = j(v + sv ). (2.25) 

□ 

Now we compile some useful properties of the maps $(t>, A), if)(v, A) and 7(1;) defined 
in the (12.151 ), (12.211 ) and (12. 141 ). These results, along with their proofs, are taken from 
ttH. 

Proposition 2.6. Let the assumptions of Theorem \2.5\ hold and let the operators $(t> , A), 
if)(v, A) and 7(f) &e defined as in (|2.15l) . (12.211) and (12.141) and to Ao and xq = v + u>o 
Z?e a* in the previous discussion. Then 

D v $(v , A ) = 0, D v ifi(v 0} A ) = 0, and ^7(^0) = 0, (2.26) 

and eac/z of these operators has the same order of differentiability as F(x, A). 
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Proof. The fact that $(t>, A), ip(v, A) and j(v) all have the same order of differentiability 
as F(x, A) follows from the definition of $(t> , A) and the Implicit Function Theorem 17 .51 
By differentiating (/ — Q)F{y + ^(v, A), A) = with respect to v we obtain 

(/ - Q)D x F(v + i/>(v, A), X)(I Xl + D^(v, A)) = 0, (2.27) 

where I Xl denotes the identity on X\ = ker(D x F(x , A )). By evaluating at (v , A ), 
where x = v + w , we obtain 

(/ - Q)D x F(x , X )D v iP{v , A ) = 0. (2.28) 

Given that D v i/j(v , A ) maps onto X 2 and (/ — Q)D x F(x , A ) is an invertible operator 
from X 2 to Yi, we have that D v ijj(x , A ) = 0. 

Then if we differentiate $(t>, A) = QF(v + ^(v, A), A) = with respect to v and 
evaluate at (t> , Ao), we obtain 

D v <f>(v , A ) = QD x F(x , Xo)I Xl = 0. (2.29) 

By differentiating (12.221) with respect to v and utilizing (12.291 ), we have 

D x $(v Q ,\ )D v j(v ) =0. 

The assumption that D\$(vo, Ao) ^ implies that 

AM^o) = 0. (2.30) 

□ 

Once we've obtained a unique solution curve (x(s), A(s)) through (x , A ), we analyze 
A(0) (where = 4-) to determine additional information about the solution curve. In 
particular, we can determine whether or not a saddle node bifurcation or fold occurs 
at (xo, Ao). This type of bifurcation occurs when the solution curve {x(s), A(s)} has 
a turning point at (a;o, Ao). The next proposition, taken from [fT3l . provides us with a 
method to determine information about A(0). 

Proposition 2.7. Let the assumptions of Theorem 12.51 be in effect. Additionally assume 
that ker(DxF(xo, Ao)) is spanned by vq. Then 



— 2 F(x(s),X(s)) 

s=0 

-.2 



(2.32) 



(2.31) 

s=0 

D x F{x , Aq)x(O) + D\F(xq, Aq)A(O) = D x F(x , X )v = 
ds 2 

D 2 xx F{x , Ao) [vq, vo] + D x F(x , A )x(0) + D x F(x , A ) A(0) = 0. 

In particular, an application of the projection operator Q defined in (|2.12l) to (12.321) 
yields 

QD 2 xx F(x , Ao) [v , vo) + QD x F(x , A ) A(0) = 0. (2.33) 
This implies that ifD\F(xo, Ao) ^ R(D x F(xq, Ao)) and 

D 2 xx F(x ,Xo)[va,v Q } £ R(D x F(x , A )), 

then A(0) ^ 0. 
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Proof. Let {x(s), A(s)} be the solution curves for F(x, A) = defined by (12.241) and 
(12.25b . Differentiating these curves we obtain 

= v + D v ip(vo, X )vo + D x ^(v 0l Xo)D v ^(v )v = v , (2.34) 

s=0 



Ts x[s ' 



^F(x(s),X(s)) 

-.2 



(2.36) 



= D v y{vo)v = 0, (2.35) 

s=0 

where the above expressions simplify as a result of Proposition 12.61 Differentiating the 
expression F(x(s),X(s)) = twice and again using Proposition 12.61 to simplify, we 
obtain 

ds 2 

D 2 xx F{x , A ) [vq, M + D x F{x , A )x(0) + D x F(x , A ) A(0) = 0, 
where 

A(0) = D^ v j(v )[vo,Vo] and £(0) = D 2 vv ^j{v , X )[v ,v ], 
by differentiating (12.341 ) and ( 12.351 ) once more with respect to s. Applying the pro- 
jection operator Q to (12.361) yields (12.331) . Then the assumptions that D\F(x , A ) ^ 
R(D x F(x , A )) and D 2 xx F(x 0l A ) [vq, M t R(D x F(x , A )) imply that A(0) ^0. □ 

The significance of Proposition 12.71 is that it gives explicit conditions that allow us to 
determine whether or not A(0) is nonzero. Heuristically, the fact that A(0) ^ means 
that A(s) has a turning point at s = 0. This means that the graph of {x(s),A(s)} 
looks like a parabola and that a saddle node bifurcation occurs at s = (cf. lfT3l ). 
If we assume that F(x, X) is at least 3-times differentiable we may expand the oper- 
ators ip(vo + sf)o,7(fo + svo)) and 7(t> + svo) about s = as a second order Taylor 
series and use (12.241) and (|2.25l) to obtain second order representations of our solutions 
{x(s), A(s)}. This is the solution approach we take to prove non-uniqueness in both the 
CMC and non-CMC cases. 

3. Main Results 

The main results of this article pertain to the following one parameter family of prob- 
lems 

-A4>+a R 4> + X 2 a T 4> 5 - a w </r 7 - 27rpe~ A 5 = 0, (3.1) 
Lw + Xb a T (j) 6 = 0. 

Here we assume that g ab is a given SPD metric with no conformal killing fields that has 
constant, positive scalar curvature. The expressions D a and A denote the derivative and 
the Laplace-Beltrami operator associated with g a b and 

Lw = -D b (Cw) ab , 

denotes the divergence of the conformal killing operator associated with g ab . Finally, we 
define 

1 1 
op = -R, a T = — t 2 , (3.2) 

H 8 12 

a w = ha + Cw) ab (a + Cw) a \ b T = \b*t. 

o o 

In general, we assume that r E C7 1,a (A^), however when we prove our CMC results we 

will additionally require that r be constant. For the remainder of this paper we assume 
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that R is a positive constant and that \a\ = (cr ab a ab ) 2 is also a nonzero constant. Notice 
that (|3.1I) has the form of (11.61) with initial data depending on A where 

r A = At, p A = e~ x p and j A = 0. 

We show that in both the CMC and non-CMC cases that solutions to (13.21 ) are non- 
unique. Our method for doing this is to apply the bifurcation theory outlined in Sec- 
tion [231 The first step in doing this is to formulate (13.11) in a way that allows us to utilize 
the framework outlined in Section! 



F((0,w),A) 



(3.3) 



3.1. Problem Setup. We now formulate (13.11) so that we can apply the Liapunov-Schmidt 
reduction. Define F((cj), w), A) by 

-A<f) + a R <j) + \ 2 a T <f) 5 - a w -7 - 27rpe~ A 5 
Lw + Xb a T <p 6 

and in the event that r is constant, define 

G(4>, A) = -A0 + a R( p + A 2 a r 5 - ^<x 2 0~ 7 - 27ipe~ x ^. (3.4) 

8 

If F(((f), w), A) = (resp. G((j>, A) = 0) for a given A, then ((0, w), A) (resp. (0, A)) 
solves Eq. d3j} (resp. Eq. (1341) ). 

We view (13.31) and (13.41) as nonlinear operators between the Banach spaces 

F((0, w), A) : C k ' a {M) © C k > a (TM) xl4 C k ~ 2 ' a {M) © C k ~ 2 ' a {TM), 

G(</>, A) : C k ' a {M) x R -)■ C fc - 2 ' a (7W). 

where k > 2. For 7^ and X = (0, w), the first order Frechet derivatives D^Gfy, A), 
D\G{4>, A), D x F(((f>, w), A) and D X F((4>, w), A) all exist. In fact, both F and G are 
A;-differentiable for any k £ N provided that 0^0. See the Appendix 17.11 for more 
information regarding Frechet derivatives. 

Now we are ready to state the main results of this paper. The first two results state 
that there is a critical density p = p c such that there exists a constant C where the 
linearizations Z)^G(0 c , 0) and DxF(((f) c , 0), 0) have a kernel of dimension one. This 
provides the basis for our final two main results where we determine explicit solution 
curves {0(s), A(s)} and {(0(s), w(s)), A(s)} to obtain our non-uniqueness results. 

3.2. Existence of p c such that dim ker(D x F((<f) c , 0, 0)) = 1. The two results in this 
section pertain to the existence of a critical energy density p = p c at which the lineariza- 
tions of the operators F and G develop a one-dimensional kernel. These results allow us 
to apply the Liapunov-Schmidt reduction outlined in Section [2~4l to analyze solutions in 
a neighborhood of ((0 C , 0), 0) and (0 C , 0). We present the theorems here without proof 
and postpone them until Section |4~2l 

Theorem 3.1 (CMC). Let D^G((j), A) denote the Frechet derivative of (13.41) with respect 
to (p. Then there exists a critical value of p = p c and a constant <p c such that when p = p c , 
Eq. (13.41 ) has a solution if and only if X > 0. Furthermore, dimker(D^G((j) c ,0))) = 1 
and it is spanned by the constant function 0=1. Moreover, we can determine the explicit 
values of p c and C , which are 

3 1 

" c = 5i^H """ ^ = (^)*- (3 ' 5) 

Proof. We present the proof in Section 14.21 □ 
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Theorem 3.2 (non-CMC). Let D X F (((/>, w), A) denote the Frechet derivative ofEq. (13.31) 
with respect to X = (0, w) and let p c and C be as in Theorem \3.1\ Then when p = p c , 

dim ker(DxF(((p c , 0), 0))) = 1 and it is spanned by the constant vector Q . 

Proof. We present the proof in Section 14.21 □ 

3.3. Non-unique Solutions to F(((f>, w), A) = when p = p c . The two Theorems in 
this section pertain to the non-uniqueness of solutions to the nonlinear problems (13.31) 
and (13.41) . Theorem l3.3l provides the explicit form of solutions to (13.41) in a neighborhood 
of the point (0 C , 0) in the CMC case. The form of this solution curve implies that a saddle 
node bifurcation occurs at (0 C , 0) and that solutions are non-unique in a neighborhood 
of this point. Theorem [3]4] provides analogous results in the non-CMC case for the point 

((0c,O),O). 

Theorem 3.3 (CMC). Suppose that r is constant. Then (13.31 ) reduces to the scalar prob- 
lem 

-A0 + a R( p + (A V - 2npe- x )4> 5 - -cx 2 0~ 7 = 0. (3.6) 

8 

When p = pc, with p c as in Theorem 13.21 then there exists a neighborhood of (<f) c , 0) 
such that all solutions to (13.61) in this neighborhood lie on a smooth solution curve 
{0(s), A(s)} that has the form 

0(s) = c + s + O(s 2 ), (3.7) 

\( s ) = h(0)s 2 + O(s 3 ), (A(0)^0). (3.8) 

In particular, there exists a 6 > such that for all < A < 5 there exist at least two 
distinct solutions 0i a ^ 02, a to (13.61) . 

Proof. We postpone the proof until Section [5TT1 □ 

Theorem 3.4 (non-CMC). Suppose r G C l ' a (M) is non-constant and let F(((f), w), A) 
be defined as in (13.31) . Then if p c and C are defined as in Theorem \3.1\ and p = p c , there 
exists a neighborhood of ((0 C , w), 0) such that all solutions to F(((j), w), A) = in this 
neighborhood lie on a smooth curve of the form 

4>(s) = C + s + ^X(0)u(x)s 2 + 0(s 3 ), (3.9) 

w(s) = ^A(0)v(x)s 2 + O(s 3 ), 

\(s) = h(0)s 2 + O(s 3 ), (A(0)^0), 

where u(x) G C 2,a (M), v(x) G C 2,a {TM) andv(x) ^ 0. In particular, there exists 
a 5 > such that for cdl < A < 5 there exist elements (0i,a? w i,a), (02 a, w 2,a) G 
C 2 > a {M) © C 2 > a (TM) such that 

F({(f>i,\, w iiA ), A) = 0, fori G {1, 2}, and (0i, A , Wi jA ) ^ (0 2 ,a, w 2 ,a)- 



Proof. We present the proof in Section [5^21 



□ 
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4. Some Key Technical Results 

4. 1 . Existence of a Critical Value p c . In this section we lay the foundation for proving 
Theorems 13.11 and |3T2| As in lUTl , we seek a critical density p c where our elliptic problem 
goes from having positive solutions to having no positive solutions. In particular, what 
we seek is a value p c such that when A = 0, then (13.31) will have no solution for p > p c 
and will have a solution for p < p c . 

When A = 0, the assumption that g a b admits no conformal killing fields implies that 



F((0,w),O) = F((0,O),O) 
Define 



-A0 + a R <p - ^<p- 7 - 2vrp0 5 = 



8 

w = 



(4.1) 



q(x) = dRX - \a 2 X 7 ~ ZnpcX 5 , (4-2) 
o 

where p c is a constant to be determined. The objective will be to determine p c so that 
q(x) has a single, positive, multiple root and then use the maximum principle discussed 
in Appendix 17.61 to conclude that if p > p c , then (14.11) will have no solution. This leads 
us to the following proposition. 

Proposition 4.1. Let q(x) be defined as in (14.21) . Then there exists constants p c > and 
4> c > such that q(x) < Ofor all x > and the only positive wot ofq(x) is 4> c - 

Proof. To determine p c , we observe that because and a 2 are constants, we simply 
need to analyze the roots of (|4.2I) as p c varies. We seek p c such that q{x) has a single, 
positive, multiple root. We observe that q(x) = if and only if 

P(X) = a RX 8 ~ ~ ^PcX 12 = 0. 

Furthermore, it is clear that each pair of roots {— xo, Xo} of the even polynomial p(x) is 
in direct correspondence with each positive root of ^(7) = a#7 2 — |o" 2 — 27rp c 7 3 , where 
7 = x 4 . Therefore, we simply need to choose p c such fhatp(7) has a single positive root. 
To accomplish this, we find the lone, local maximum of ^(7) and require it to be a root 
of ^(7). We have that 



and 



= p'hf) = 2a#7 — 67rp c 7 2 =^ 7 C = — — is a local max, 

3vrp c 



\ 1 1 / 
a R \ 1 2 o _ / a R 



°^ ( ^^H^J -r 2npc {^p-) (4 - 3) 

4-JoWpj) rI 

27tt 2 P 2 Pc 24^1^1^' 

□ 



The next result follows immediately from the previous analysis but will be useful going 
forward. 
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Corollary 4.2. Define the constants 

Rl and c =f^V. (4-4) 



24V3|cr|7r " \3irp. 
Then if 

q(x) = a RX - g a2 X~ 7 - 2np c x 5 , 
it follows that q ((p c ) = q'(4> c ) = 0. 

Proof. This follows immediately from the proof of Proposition 14. II or by direct compu- 
tation. □ 

Now we show that p c is a critical value of (|4.1I) . 

Proposition 4.3. Let p(x) G C(M). Then the constant p c defined in Corollary \4.2\ has 



the property that Eq. (14.11) has a positive solution if < p < p c and has no positive 
solution if p > p c . 

Proof. Let q(x) be defined as in Corollary 14 .21 If <p > solves (14.11) . then 

= M> - la 2 r 7 ~ 2vrp0 5 = f{x, <P). (4.5) 

o 

We observe that if p > p c , then p = inf ^.m p > p c and for x > 0, 

/(z, X) = - T^X^ 7 - 2?rpx 5 < orX - l^X^ 7 ~ ^PX 5 < <?(x)- ( 4 -6) 
o o 

Therefore if p > p c , (14.51) and (14.61) imply that any positive solution <p to (|4.1I) satisfies 

A0 = /(x,0) <g(0) <0. 

An application of the maximum principle (|7.6I) implies that if p > p c , then (14.11) has no 
solution. 

To verify that (14.11) has a solution if p < p c , first observe that Corollary 14.21 implies 
that 

*,=r^y =(-?-y, (4.7) 

solves Eq. (14.11) when p = p c . If p < p c , the properties of imply that the polynomial 

<?i(x) = aRX - l° 2 X~ 7 - 2ttPX 5 , p = sup p(x), 

will have two positive roots xi < X2- Therefore, any (j) + satisfying < xi < 4>+ < Xi 
will be a positive super- solution to (14.11) given that 

f{x, X) > <?i(x) = aRX ~ \° 2 X~ 7 ~ 2vrpx 5 - 

Similarly, we may choose a positive sub-solution 0_ < + to (|4.1I) by choosing any 
sufficiently small 0_ satisfying < 0_ < % 3 , where % 3 is the lone positive root of 

taix) = a R x - \ a<2 X~ 7 - 

We can then apply the method of sub- and super-solutions outlined in Section 17.2.21 to 
solve (14~TT). □ 



The next result extends Proposition 14.31 to the case when A ^ and indicates that p c is 
also a critical value for the decoupled problem (13.41) . 
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Corollary 4.4. Let p(x) G C(M) and suppose that r is a constant and that 



Pc 



2AV3\a\n 

There exists an e > such that there is no positive solution to (13.41) if p > p c and 
— e < A < 0, and there exists a positive solution to (13.41) if0<p<p c and < A < e. 
Finally, if p = p c and A is sufficiently small, then (13.41) has a solution if and only if\>0. 

Proof. Again, we observe that if cj) > solves (13.41) , then 

A<f) = a R <p + A 2 a T 5 - \a 2 <p~ 7 - 2vrpe-> 5 = f(x, 0, A). (4.8) 

8 

Let q(x) be as in Corollary 14.21 and define 

pi(x, A) = a R x + AVx 5 - \<y 2 X^ - 2vrpe _A x 5 , 

o 

where p = mf xeM p(x). It is clear that f(x, (f), A) < pi((f), A) for any 4> > 0, and for 
A < and p > p c we have that 

Pi (x, A) = a RX + X 2 a rX 5 - -a 2 X ' 7 ~ 271-pe" V (4.9) 

o 

< a RX + (A 2 a r - 2vrp c + 2vrp c A + o(\ 2 )) X 5 - \^X^ 
= q(x) + (AV + 2vrp c A + o(A 2 ))x 5 = q(x) + #(A)x 5 . 

Here we observe that g(\) — > as A — > 0, and for |A| sufficiently small, g(\) < if 
A < 0. By Proposition 14.11 we know that if X > then q(x) < 0. So Eq. (|4.9I) implies 
that if p > p c and A < is sufficiently small, then f(x, X , A) < Pi(x> A) < 0, and the 
maximum principle then implies that (13.41) will have no solution. 
If P < Pc, then define 

p 2 ( X , A) = a R x + A 2 a r X 5 - \ a2 ^ 7 ~ 27[ P e ~ X X 5 , 

where p = sup^.^ p(x). It is clear that f(x, X -> A) > P2(x A) for all X > 0, and for 
A < we have 

MX, A) = a R X + AVx 5 - \o 2 X~ 7 ~ ^P^\ h (4.10) 

o 

> a RX + (A 2 a r - 2vrp c + 2vrp c A + o(A 2 ))x 5 - -<J 2 X~ 7 

o 

= q( X ) + (X 2 a T + 27rp c A + o(A 2 )) X 5 = q( X ) + 9Wx"- 

Again, g(X) — > as A — > and g(\) > for A > sufficiently small. Therefore if 
X > 0, Eq. (14.101) implies that f(x, x, A) > P2<X A) > <?(x) if A > 0. The properties 
of q{x) specified in Proposition 14. 1 1 imply that that for any A > 0, either P2(x A) has 
a single positive root Xo and P2(x A) > for all x > Xo> or p 2 (x A) has two distinct 
positive roots. This implies that if A > we can find a positive super-solution <p + to (I3.4I ). 
If A = we take 0+ = <p c to be a super-solution where C is defined in Corollary 14.21 
Similarly, we can also find a positive sub-solution 0_ satisfying 0_ < + by choosing 
any sufficiently small < 0_ < xo, where xo is the unique positive root of 

r(x, A) = a R x + \ 2 a T x 5 - \° 2 X~' '• 
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The method of sub-and super-solutions outlined in Section 17.2.21 then implies that if 
p < p c and A > 0, then (13.41 ) has a solution. 

Finally, we observe that if p = p c , then we have that 

f{x, x ,X) = q{x) + gWx 5 , 

where / and g are the same as above. Therefore, when A is small and p = p c , we 
can apply the above analysis to conclude that (|3.4[) will have a solution if and only if 

A > 0. □ 

Remark 4.5. We note that the negative sign in front of the term 2irpx 5 in the polynomial 

I 

?(X) = a RX ~ g^X" ~ 27 W\ 

played an essential role in allowing us to determine our critical density p c and critical 
solution (p c . If this term were positive, then q(x) would be monotonic increasing for 
X > 0, and we would not be able to find a positive (p c and p c so that q((f> c ) = and 
q'(<p c ) — 0. As we saw in Corollary \4.4\ and Proposition \4.3\ these properties of q(x) 
played an important role in the existence of solutions to Eq. (13.41) and Eq. (|4.1I) . Later in 
this article, we will also see that these properties ofq(x) play an important role in our 
non-uniqueness analysis by allowing for the kernel of the linearization of F (((f), w), A) 
and G((f), A) to be one-dimensional. These facts further emphasize the role that terms 
with the "wrong sign" (cf. IfTTl ) have in the non-uniqueness phenomena associated with 
the CTS, CTT and XCTS formulations of the Einstein constraint equations. 

4.2. Existence of a One Dimensional kernel of D X F((4> C , 0), 0) when p = p c . In the 

previous section we proved the existence of a critical density p c that affected whether 
Eq. (14.11) and Eq. (13.41) had positive solutions. We now show that when p = p c , the 
linearization of both (13.41) and (13.31) develops a one-dimensional kernel. 

We first calculate the Frechet derivatives DxF(((f>, w), A) and D^Gfy, A). To com- 
pute these derivatives, we need only compute the Gateaux derivatives given that the G- 
derivatives are continuous in a neighborhood of (((f) c , 0, 0). See [20] and Remark [7721 
Therefore, 

D x F(((f> c , 0), 0) = j f F((<f>c + t(f>, tw), 0) 

where O, w) e C 2 ' a (M) © C 2 ' a (T M) satisfies \\((f>, w)\\c^{M)®c^{TM) = !• 
So for a given ((0, w), A), the Frechet derivative 

D x F((<p, w), A) : C 2 ' a (M) © C 2 ' a (TM) ->• C°' a (M) © C°' a (TM), 

is a block matrix of operators where the first column consists of derivatives of 

F(((f>, w), A) with respect to <\> and the second column consists of derivatives with respect 

to w. This implies that 



D x F(((f>,w),\) 
where 



-A + a R + 5A 2 a r 4 + 7a w </r 8 - 10np c e~ x (p 4 L 
6AK(/> 5 L 



(4.11) 



hh = h((f>,w)h = ((Cw) ab (Ch) ab + a ab (Ch) ab ) , (4.12) 

and C is the conformal Killing operator. Similarly, in the CMC case the map 

Z^G(0, A) : C 2 > a (M) -> C°> a (M), 
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has the form 

D^G((p, A) = -A + a R + 5A 2 a T 4 + ^a 2 ^ - 10np c e- X (p 4 . 

8 

We now make some key observations about (14.111) . 

Proposition 4.6. Let <p c be as in Corollary \4.2\ Then F((0 C , 0), 0) = and 
DxF((4> c , 0), 0) has the form 



(4.13) 



D X F(0 C ,O,O) 



-A L 
L 



(4.14) 



where L : C k > a {TM) -»■ C*"^"^) is defined by 

1 



06 



L(0 C) 0)/i = L/i = ~C 7 ^(^) 
anJ £ is the conformal killing operator. 

Proof. By Corollary 14.21 it follows that <p c is a root of the polynomial 



and also a root of 



q{x) = a RX - 7 ~ 2np B )£ 



q'ix) = a R + -a^ 8 - lOnp cX 4 . 



(4.15) 



This implies that F((0 C , 0), 0) = and that Eq. (I47TTT) reduces to (1041) when 

((0,w),A) = ((0 c ,O),O). 

Remark 4.7. Corollary \4. 21 implies that (|4.13l) reduces to 

D+G{<f> e ,0) = -A, 

in the CMC case. Therefore dimker{D ( f l G{(j) Cl Q))) = 1 and it is spanned by the constant 
function = 1. 

Corollary 4.8. Letting Hi = L 2 (M) and U 2 = L 2 (T M), the Hi © U 2 -adjoint of 
D x F(((f) c , 0), 0) has the form 



□ 



(4.16) 



(D x F(4> c ,0,0)y 



-A 

£ L 



(4.17) 



(4.18) 



where L : C k ' a (M) C k ~ 1 ' a (T M) is defined by 

tu = D\-^- c 7 ua ab ). 

Proof. Let (ui, Vi) and (« 2 , v 2 ) both be elements of C 2 (M)@C 2 {T M). Then given that 
both -A and L = -D b (C) ab are self-adjoint with respect to the L 2 (M) and L 2 (T M) 
inner products, it follows that 



D x F(((j) c ,0),0) 



Ui 




U 2 


/ J M 


. V1 . 


5 


. V2 . 



{—U\Au 2 + vi • Lv 2 + ]Lv 1 u 2 )dV g: 

(4.19) 
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where dV g is the volume element associated with g ab and Lvi = — \4>~ 7 o ab (C\i) ab . 
Given that the negative divergence of a (0, 2) tensor and the conformal killing operator 
C are formal adjoints (see 112010 . we have that 



M 



M 



M 



1 



^ 2 c -V a6 (£ Vl ) a6 ) dV g 



(4.20) 



1 



D\-u 2 <p-'a ab )-w l )dV g 



hu 2 • VidVg 



M 



Ml 




u 2 




j 


. V2 . 



Therefore, 

D x F{(<j> c ,0),0) 

{—u\Au 2 + vi • Lv 2 + Lw 2 • \i)dV g 

M 



Corollary 4.9. D x F(((f) c , 0, 0) has a kernel of dimension 1 that is spanned by 
(D x F((f) c , 0, 0))* also has a kernel of dimension one that is spanned by 



(4.21) 



Ml 




' -A 


" 




u 2 




1 


L 


L 







□ 

and 



Proof. We solve for 



G C 2 ' a (M) © C 2 > a (TM) such that 



u 




—A L " 




u 




" " 


V 




L 




V 








£>*F((0 C ,O),O) 



Given that g ab admits no conformal killing fields, we must have that v = 0. This implies 
that 

= —Au — -<f)~ 7 (o- a b(£v) ab ) = —Au ==>• u is a constant. 



Therefore 



spans ker(£> x F((0 c , 0),0). 



Similarly, we solve for 



such that 



u 




" -A 


" 




u 




" " 


V 




L 


L 




V 








OD x F((0 c ,o),o)y 

This implies that u is a constant and that 

= Lu + Lv = V b (^(f) c ucr ab ) + Lv = ^(p c uV b a ab + Lv. 
Given that a ab is divergence free, we have that V h a ab = 0, which implies that v = 0. 



Therefore 



spans ker( J D x F((0 c ,O),O)*) 



□ 



We can now prove Theorems 13.11 and 13 .21 The proofs are an immediate consequence 
of the preceding results, but we summarize them here in the proof for convenience. 

4.3. Proofs of Theorems 13.11 and I3.2t Critical Parameter and Kernel Dimension. 

Proposition 14.3 1 implies the existence of critical values 



Pc 



(— 



and <p c 



an 
3np c 
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such that if 



q{x) = a RX 



■ 2 -7 

r 0- x 



then q((f) c ) = q'((f) c ) = 
CMC case reduces to 



0. By Remark 14771 we have that the linearization (14.131) in the 
-A. This proves Theorem 13. II Similarly, in Proposition l4T6l we 



explicitly determined D x F(((j) c , 0), Cf 
kernel spanned by the constant vector 



and in Corollary 14.21 we showed that it has a 



This proves Theorem 13 .21 



4.4. Fredholm properties of the operators D x F(((f> c , 0), 0) and D^G^cO). Now 
that we have shown that the linearizations D x F(((j) c , 0), 0) and D^G((j) c , 0) have one- 
dimensional kernels, we are almost ready to apply the Liapunov-Schmidt reduction. Re- 
call from section[2]that a key assumption in this reduction was that the operator be a non- 
linear Fredholm operator. Therefore, to apply this reduction in the CMC and non-CMC 
cases we must show that the operators D^G^c, 0) and D X F((4> C1 0), 0) are Fredholm 
operators between the spaces on which they are defined. In particular, we need to show 
that D^G^c, 0) is a Fredholm operator between the spaces C 2 ' a (Ai) and C 0,a (Ai) and 
that the operator D x F(((f) c , 0), 0) is a Fredholm operator between C 2 ' a (M)®C 2 ' a (T M) 
and C°' a (M) ®C°' a (TM). 

In the CMC case, we have that D < j ) G{4> c , 0) = — A. It is well known that this operator is 
a Fredholm operator between the Hilbert spaces H 2 (A4) and L 2 (A4) [TTOl . Furthermore, 
—A is a Fredholm operator between the subspaces C 2 ' a (Ai) and C 0,a (M.) because of 
the regularity properties of the the Laplacian and the fact that these spaces continuously 
embed into the Hilbert spaces H 2 (A4) and L 2 (A4). See Appendix 17.2.31 for a more 
detailed discussion of these facts. 

Letting L = —A, we regard L = L* as operators from H 2 (M) — > L 2 (M). The 
Fredholm properties of these operators allow us to make the following decompositions 
that are orthogonal with respect to the L 2 -inner product: 

L 2 (M) = R(L*) © ker(L) (4.22) 
L 2 (M) =i?(L)©ker(L*). 

In this case, these decompositions are the same given that L is self-adjoint. Therefore if 
we regard C 2 ' a (M) and C 0,a (M) as subspaces of L 2 (M), then we may use (14.221 ) to 
obtain the following decompositions 

C 2 ' a (M) = (R{L*) n C 2 ' a {M)) © ker(L), (4.23) 
C Q ' a (M) = {R{L)nC°' a {M))®kev{L*), 

which are also orthogonal with respect to the L 2 -inner product. See Appendix 17 . 2 . 3 1 for 
further details. 

It is not as clear that the operator D x F(((p c , 0), 0) is a Fredholm operator between the 
spaces C 2 > a {M)®C 2 > a {TM) and C°' a (M ) © C°' a (T M ) . For the sake of completeness, 
we briefly discuss this point. As in Appendix 17.2.31 we first show that D x F((p c , 0), 0) 
is a Fredholm operator from the Hilbert space L 2 (Ai) © L 2 (TAi) to itself, where we 
consider the domain of definition of -Dx-F((0 c , 0), 0) to be H 2 (M) © H 2 (T M). Indeed, 
the operator D x F(((j) c , 0), 0) induces the bilinear form 

B((u u v 1 ),(u 2 ,v 2 )) : (H\M)@H l (TM)) x (H l (M) © ^(TM)) -> R, 
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where (■, ■) is the inner product associated with L 2 (Ai) © L 2 (TM) and 
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B((ui, Vi), (it 2 , v 2 )) 



-A L 
L 





Ml 




«2 




Vl 


5 


v 2 



(4.24) 



Paralleling the discussion in 17.2.31 we first show there exists constants C, c > such that 

B((u, v), (u, v)) + c((ix, v), (u, v)) > C|| (u, v) \\ 2 H i {M)(BH i {TM y 
Let c > be a constant to be determined. Then 

B{{u, v), (it, v)) + c<(«, v), (u, v)) (4.25) 
= y (D a uD a u - ^0 C -V a6 (£^) 06 + (£i;) a6 (£i;) afe + cu 2 + cv a v)j dV g 



> 



M 



D a uD„u 



1 



16ce 



-u — e< 



C (a ab (Cv 



,ab\2 



(Cv) ab (Cv) ab + cu 2 + c^ a ) rfV £ 



' 9' 



where the above inequality follows from an application of Young's inequality. The 
Schwartz inequality and the definition of £ then imply that 

< 7 a6 (/:i;) a6 = ( < 7,/:v) ff <CH|Dv|. 

Therefore 



-14, 



,ab\2 



<c(e) 



11,2' 



M 



where c(e) -> as e -> 0. Combining (14.251) and (14.261) we have that 

B((u, v), (it, v)) + c((it, v), (u, v)) > 



(4.26) 



(4.27) 



'1 - cfel 



|v||?, 2 + \\Du\\l,2 + (c - y^)II m IIo,2 > C(||v||? i2 + Hk) 



where the final inequality holds by choosing e sufficiently small and c sufficiently large. 
The above discussion tells us that the bilinear form 

B((u, v), (u, v)) + c((it, v), (u, v)) 

is coercive on i/ 1 (A1) © H l {TM). The Lax-Milgram theorem implies that the problem 



it 






V 







( J D x F((0 c ,O),O) + c/) 

has a unique weak solution (it, v) 6 iJ x (A^) © H l {TM) for each (/, g) 6 
L 2 (.M) © -^ 2 (T A^), and elliptic regularity gives us that (it, v) e H 2 (M) © H 2 (T M). 
Therefore we conclude that the operator DxF((cj) c , 0), 0) + cl is a bijection between 
H 2 (M) © H 2 (T M) and L 2 (M) © L 2 (T M). We are the able to conclude that 

(D x F(((j) c , 0), 0) + ciy 1 exists and is compact. 

Paralleling the discussion in Appendix 17.2.31 we can then conclude that the operator 
D x F(((f) c , 0), 0) is a Fredholm operator between H 2 (M) © H 2 (T M) and 
L 2 (M) © L 2 (TM). Using the fact that l7°> q (.M) © C°' a (TM) embeds continuously 
into L 2 (.M) © L 2 {J~M) and invoking classical Schauder estimates, an argument similar 
to the argument in 17.2.31 implies that DxF(((f) c , 0), 0) is Fredholm operator between 
the spaces C 2 ' a (M) © C 2 ' a (TM) and C°' a (M) © C°' a (T.M). By applying the same 
argument to DxF(((p c , 0), 0)*, we can also conclude that this operator is a Fredholm 
operator between C 2 > Q (.M) © C 2 ' a (TM) and C7°> Q (.M) © C7°< a (T.M). 

If L = D x F(((j) c , 0), 0), then the fact that both L, L* are Fredholm operators from 
H 2 {M)®H 2 {TM) -> L 2 {M)®L 2 {TM) allows us to decompose L 2 {M ) ®L 2 {TM) 
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as in ( 1522b . Therefore, regarding C 2 ' a (M) © C 2 ' a (TM) and C°' a (M) © C°> a (TM) 
as subspaces of L 2 (.M) © L 2 (TA4), we obtain the following decompositions that are 
orthogonal with respect to the L 2 (M) © L 2 (TM.)- inner product: 

C 2 ' a (M) © C 2 ' a (TM) = ker(L) © (R(L*) n (C 2 ' Q (.M) © C 2 '°(T.M))), (4.28) 

C°' Q (7W) © C°' Q (T.M) = ker(L*) © (R(L) n (C°' Q (.M) © C°' Q (T.M))). 

In the above decomposition, L = D x F(((f) c , 0), 0) and ker(L), R(L), ker(L*) and 
i?(L*) are all regarded as subspaces of L 2 (.M) © L 2 {TM). 

5. Proofs of the Main Results 

5.1. Proof of Theorem I3.3t Bifurcation and non-uniqueness in the CMC case. We 

are now ready to prove Theorem 13.31 In the CMC case, our system (13.11) with p = p c 
reduces to 

G((f>, A) = -A0 + a R <p + \ 2 a T <f - -a 2 4>~ 7 - lTxp c e~ x ^. (5.1) 

8 

To prove that solutions to (|5.1I) are non-unique, we will apply the Liapunov-Schmidt 
reduction outlined in Section 12741 and then invoke Theorem [23] and Proposition [277] 

By Theorem 13. II and Remark 14771 we know that D^G^cO) = —A. It follows that 
dim ker(_D^G(0 c , 0)) = dim ker(D^G((p c , 0)*) = 1, where both spaces are spanned by 
0=1. 

Using the notation from Section I2.4L we can apply the Liapunov-Schmidt Reduction, 
where vq = 1 is a basis of ker^D^G^c, 0)) = ker(D^G(0 c , 0)*). By the discussion 
in Section 1474] and appendix 17.2.31 we can decompose X = C 2,a (Ai) = X\ © X 2 and 
Y = C^ a (M) =Y X @ Y 2 , where 

X x = ker(r^G(0 c , 0)), X 2 = R{D^G(<p c , 0)*) n C 2 > a (M), (5.2) 

Y 1 = R(D (/ ,G(<f> e ,Q))nC ' a (M), and Y 2 = ker(D^(0 c , 0)*). 

Letting P : X — )■ X x and Q : Y — > Y 2 be projection operators as in Section 12.41 and 
writing <p = P<t> + (I — P)<f> = v + w, the Implicit Function Theorem applied to 

(I -Q)G(v + w,X) = 0, (5.3) 

implies that w = ip(v,X) in a neighborhood of (0 C , 0) and = ip((j) c ,0). Plugging 
, A) into 

QG(v + w,X) = 0, 

we obtain 

$(v,\) = QG(v + ij(v,X),\) = 0. (5.4) 

All solutions to A) = in a neighborhood of (0 C , 0) must satisfy Eq. (15.41) . 
We now observe that D\G(4> C , 0) = 2 r xp c <\)\ ^ 0. This implies that 

£>A$(0 C , 0) = QP A G(0 c , 0) = 2vrp c 0^ ^ 0, (5.5) 

given that Q is the projection onto Y 2 and Y 2 is spanned by the constant function 1. 
The Implicit Function Theorem applied to Eq. (|5.4I) implies that there exists a function 
7 : U± — > V\ such that U\ C Xi, Vi C R and 7(f) = A in a neighborhood of C with 

7(&) = 0. 

Therefore (15.41) becomes 

g(v) = QG(v + i>(v n (v)) n (v)), (5.6) 



NON-UNIQUENESS OF SOLUTIONS TO THE CONFORMAL FORMULATION 23 

and by writing v = s + <j> c , which we can do for s £ (—5, 5) with 5 > sufficiently 
small, we obtain 

g(s) = QG(s + C + i)(s + <f> , 7 (s + C )), 7(s + C )) = 0. (5.7) 

This implies that solutions to G(4>, A) = are given by g(s) = in a neighborhood of 
(4> c , 0), where 

4>{s) = s + c + ^(s + c , 7 (s + c )), (5.8) 
A(s) = 7 (s + c ) 

determine a differentiable solution curve through (0 C , 0). 

Equation (15.81 ) gives us a fairly explicit representation of the continuously differ- 
entiable curve {0(s),A(s)} provided by Theorem 12.51 However, by applying Propo- 
sition 12.71 we can determine that A(0) ^ to obtain even more information about 
{4>(s), A(s)}. We observe that 

D 2 H G{cj> c , 0) [v , vo] = -7a 2 <p- 9 - AOirp^l + 0. (5.9) 

Therefore 

-7a 2 c " 9 - 40vrp c ^ e Y 2 D^G(&, 0)[t) , «o] £ R{D„,G{(j> c , 0)) = F x , 

given that Yi _L Y~ 2 . Proposition 12.71 implies that A(0) ^ and that a saddle node 
bifurcation occurs at (0 C , 0). 

We now combine (15.81 ) and the fact that A(0) ^ to obtain a more explicit representa- 
tion to the solution curve {(j>(s), A(s)} in a neighborhood of (0 C , 0). Define the function 

/( S ) = ^( S + C , 7 ( S + C )). (5.10) 
Then by Propositions 12.61 and 12.71 we have that 

/(0) = 0, and A(0) = 7 (0 C ) = 0, (5.11) 



A(0) = ±\{s) 
as 



= D vl (<f) c ) = 0, 

s=0 



= D v ^(cj) c , 0) + £>aV>(0c 0) A/y(<k) = 0. 

Therefore the function f(s) = 0(s 2 ). By computing a Taylor expansion of A(s) about 
s = and using Eq. (15.111) and Eq. (15.81) . we find that for s £ (—5, 5), 

<f)(s) =<j) c + s + 0(s 2 ), (5.12) 
A(,) = iA(0), 2 + O( S 3 ), 

where A(0) ^ 0. 

Based on the form of 0(s) and A(s) in Eq. (|5.12l) . there exists a <5' £ (0, 5) such that 

<f>(s) < 0, A(s) > for all s E [S',0), and > 0, A(s) > for all s E (0,5'}. 
Letting M = min{Mi, M 2 }, where 

Mi = sup A(s) and M 2 = sup A(s), 
se[-(5',o] se[o,<5'] 

the Intermediate Value Theorem then implies that for all Ao £ (0,M), there exists 
Si,s 2 E [—5', 5'], si 7^ s 2 , such that A(si) = A(s 2 ) = A . Based on how we chose 
5', we also have that <f>(si) ^ 0(s 2 ). This completes the proof of Theorem l3.3[ 
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5.2. Proof of Theorem |3.4t Bifurcation and non-uniqueness in the non-CMC case. 

In this section we will show that solutions to F((4>, w), 0) = for the full system 



F((0,w),A) 



— A0 + an4> + \ 2 a T (, 



k-7 



Lw + Xb 



2irpe 



(5.13) 



are non-unique, where r 6 C 1,a (Ai) is a non-constant function. Our approach is similar 
to that of the CMC case: we apply a Liapunov-Schmidt reduction to Eq. (|5.13l) to deter- 
mine an explicit solution curve through the point (((f> c , 0), 0). The form of this curve will 
imply that solutions to the system (15.131) are non-unique. 

By Proposition 14.61 we know that kerD x F(((j) c , 0), 0) takes the form 



D x F{<p c , 0,0) 



-A L 
L 



where hh 



4^c 



7 a ab {Ch) ab . Corollary gives us that ker(D x F(((j) c , 0), 0)) and 
ker(D x F(((f) c , 0), 0)*) are spanned by v = * . 

Using the notation from Section 12.41 we apply the Liapunov-Schmidt Reduction. By 
the decomposition (14.281) . we have that 

X = C 2 ' a {M) © C 2 ' a (TM) = Xi 



and 
where 



X, 



Y = C^ a (M) © C°' a (TM) = Y l © Y 2 , 
X 1 = ter{D x F{{<f> c , 0),0)), 

X 2 = R(D x F((<p c , 0), 0)*) n (C 2 > a (M) © C 2 ' a (TM)), 
Y 1 = i?(D x F((0 c , 0), 0)) n (C°> a (M) © C°' a (TM)), 
y 2 =ker( J D x F((0 c ,O),O)*). 



(5.14) 
(5.15) 
(5.16) 
(5.17) 



Let P : X — > X\ and Q : Y — > Y 2 be the projection operators defined using v as in 
Section [2~4l Then by writing 



v + y, 



' " 


= P 


' " 


+ (I-P) 


' " 


W 




W 


W 



the Implicit Function Theorem applied to 

(I-Q)F(v + y,X) = 0, 
implies that solutions to F (((/>, w), A) = satisfy 

$(v,\) = QF(v + ij(v,X),X) =0 
in a neighborhood of ((0 C , 0), 0), where y = ip(v, X) in this neighborhood and 

(o,o) = v((0 c ,o),o). 

We now observe that 

O-rrn A$ 

7J A F((0 C ,O),O) = 



b a T <P 6 c 



due to the fact that 



E Y 2 and Yi JL Y 2 . 



(5.18) 
(5.19) 
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This implies that 
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D A $((0 C , 0),0) = QD x F((<f> c ,O),0) 







(5.20) 



given that Q is the projection onto Y 2 . The Implicit Function Theorem again implies that 
there exists a function 7 : U% — > V\, where (0 C , 0) G XJ\ C Xl, V\ C K. and 7(11) = A in 
U\ with 7(0 C ) 0) = 0. Using this fact, Eq. (15.191 ) becomes 



and by writing 





= QF(v + ij(v, 


7(«))>7( 


v)) = 






" 1 " 






V = 


(s + <p c )v = s 





+ 






0. 



(5.21) 



for s E (—6, S) with S > sufficiently small, we then obtain 

g(s) = QF (sv + <p c v + ip(sv + (j) c v , 7(sw + c t) o )), l(sv + (p c v )) = 0. (5.22) 

This implies that solutions to F(((j), 0), A) = in a neighborhood of ((0 C , 0), 0) satisfy 
g(s) = 0, where 



4>{s) 
w(s) 



+ 







1 





+ 







.7 s 



+ 







(5.23) 



X(s) = 7 [s 



, — 1 









+ 






determine a smooth solution curve through ((0 C , 0), 0). 

As in the CMC case, we seek additional information so that we can further analyze 
the solution curve (15.231) . Now we apply Proposition [277] to determine information about 
A(0), and then we will expand the function 

/(s) = V((s + 0c)uo,7((s + ^c)«o)) (5-24) 

as a Taylor series to obtain a more explicit representation of {(<j>(s), w(s)), A(s)}. 
Taking the second derivative of F(((f>, w), A), we have that 



D 2 xx F((0 c , 0),0)[v ,v ) 



-7a 2 (f)- 9 - 40np c 4> 3 c 




eY 2 . 



(5.25) 



Given that the vector (15.251) lies in Y 2 and Y x ± Y 2 , 

D xx F((cj )c ,0),0)[v ,v ]^Y 1 . 

We can therefore apply Proposition [277] to conclude that A(0) 7^ 0. 

Our next goal is to expand the function f(s) as a Taylor series about 0. In order to do 
this, we use (|5.23l) . Propo sition 12 . 61 and the fact that A(0) 7^ to obtain information about 
coefficients in this expansion. In particular, the objective is to determine information 
about the coefficient of the second order term in the expansion of f(s). 

By differentiating 

(I-Q)F(v + 1>(v,\),\) = 0, 
with respect to A and evaluating the resulting expression at ((0 C , 0), 0), we obtain 

(/ - Q)D x F((<f) c , 0), 0)D x ^((<p c , 0), 0) + (J — Q)D x F((<f) c , 0), 0) = 0. (5.26) 

Given that 

D A F((0 C ,O),O)= \ 2 l^f 
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and Q is the projection operator onto Y 2 , which is spanned by 



1 




we have that 



(/ - Q)D\F(((f) c , 0), 0) = 
Equations (15.271) and (15.261) imply that 

(/ - Q)D x F{(<f> c , 0), O)£W((0 C , 0), 0) 







(5.27) 



(5.28) 



Given that D X F((4> C , 0), 0) has the form (14.141) and the operator L is invertible, Eq. 
(15.281) implies that 



D A V((0 C ,O),O) 



with v(a;) 7^ 0. 



(5.29) 



As we shall see, this fact implies that w(s) has quadratic terms in s. 

We have one last piece of data left to determine the coefficient of the second order 
term in the Taylor expansion of f(s). Differentiating (/ — Q)F{v + ip(v, A), A) = 
twice with respect to v, evaluating at ((0 C , 0), 0) and applying the resulting bilinear form 
to vq, we obtain 



(I-Q)D xx F((4> c ,0),0)[v ,Vo}+ 

(I - Q)D X F(((/> C , 0), O)Z^((0c 0), 0)[v Q , v ] 



(5.30) 



0. 



By Eq. (15.251) we know that D XX F((4> C , 0), 0)[v Q , v ] G F 2 -Because (7 - Q) projects 
onto Yi and Y\ _L Y 2 , we have that 



(I-Q)D 2 xx F((cj )c ,0),0)[v ,vo} = 0. 



(5.31) 



Equations (|5.31l) and (15.301) and the invertibility of (I—Q)D x F(((p c , 0), 0) as an operator 
from X 2 to Y"i imply that 



D^((0c,O),O)[uo,Uo] = O. 



(5.32) 



This was the final piece of information that we needed to to determine the second order 
expansion of f(s). 

We now expand the function f(s) in Eq. (15.241) about s = 0. We have that 



/(O) = V((0 C ,O),O) 








(5.33) 



/(0) = A,V((0 c , 0), 0)«q + D A V((0 c , 0), 0)) A/rOk, 0){) 








/(0) = L>^((0 C , 0),0)[« ,uo] +^((^,0),0)[t)o,D„7(^,0)t5o] 

+ £>LV((0c, «)> °) [A,7(&, 0)«o, t> ] + ^((0c 0), O)^7(0c 0) [t> , £0] 

+ D 2 xx ^((<l>c, 0), 0) [A, 7 (&, 0)«o, A,7(0o 0)«o] 



L> A ^((0 c ,O),O)D^7(0 c ,O)h^ o ] = ^((0 C ,O),O)A(O) ^ 
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where /(0) simplifies as a result of Proposition 12.71 Eq. (15.291) and Eq. (15.321) . which 
imply 



A,V((0 c ,o),o) = o, 

D 2 v J((<p c ,0),0)[v ,v ] 
Therefore it follows that 



A/y( 



o) = o, 



(5.34) 








^((0 C ,O),O)^ 








/(s) = -(,D A ^(0 c {) o ,7(0 c {)o))A(O)) S 2 + O( S s 



~u(z)A(0) 
_ ±v(z)A(0) 



s 2 + 0(s 3 ), (5.35) 



where we identify D x i/j(((j) c , 0), 0) with the vector in C 2 - a (.M) © C 2 ' a (TM). 

By Eq. (15.291 ) we have that v(sc) 7^ and expanding out A(s) as a second order Taylor 
series about s = we obtain 



AW 



(5.36) 



Putting together (15.231) , (15.351 ) and (15.361) we find that solutions to F((<f), w), A) = in a 
neighborhood of ((0 C , 0), 0) take the form 

0( s ) = C + s + iA(0)?i(x)s 2 + 0(s 3 ), (5.37) 



w(s) = ^A(0)v(x)s 2 + O(s 3 ), 
A( S ) = iA(0) S 2 + O( S 3 ), 



(5.38) 
(5.39) 



where s e (—5, 5) for sufficiently small 5 > 0. 

By analyzing the solution curve (|5.37I) - (I5.39I ) as we did for the curve (15.121) in the proof 
of Theorem l3.3l we can conclude that solutions to the system (15.131) are non-unique. This 
completes the proof of Theorem 13 .41 

6. Summary 

We began in Section [2]by introducing our notation for function spaces and presenting 
the basic concepts from functional analysis and bifurcation theory that we used through- 
out this paper. In particular, we gave an outline of the Liapunov-Schmidt reduction that 
was the basis of our non-uniqueness arguments. Then in Section [3] we presented our main 
results, which consisted of the existence of a critical solution where the linearizations of 
our system 

+ a R (f) + X 2 a T (f) 5 - a^^r 1 - 2npe^ x (p 5 
Lw + Xb a T <p % 

developed a one-dimensional kernel and non-uniqueness results for solutions to 
F (((/), w), 0) = in both the CMC and non-CMC cases. We then set about proving these 
results in the following sections. In Section |4~T1 we showed that in the CMC case there 
exists a critical density p c for the operator 



F((0,w),A) 



(6.1) 



G((f), A) = — A0 + clr4> + A 2 a r — a w 7 — 2-npe 



-Xj.5 



(6.2) 



This density satisfied the property that if |A| was sufficiently small, then p > p c and 
A < implied that there was no solution to G((f>, A) = 0, and if p < p c and A > 
then there was a solution. This result provided the foundation in Section [4~2l for showing 
that the linearization of (16.11) developed a one-dimensional kernel. Then in Section 14.41 
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we briefly discussed the Fredholm properties of the linearized operators D x F(((f) c , 0), 0) 
and D^G((f) c , 0) on the Banach spaces on which they are defined. 

In Section I5T1 we proved the first of our non-uniqueness results. We showed that in 
the event that the mean curvature was constant, the decoupled system (16.21) exhibited 
non-uniqueness. This was indicated by the fact that the solution curve through the point 
(0c 0) had the form 

<f>{s) =<j) c + s + 0(s 2 ), (6.3) 

X(s)= 1 -X(0)s 2 + O(s 3 ), 

which implied that a saddle-node bifurcation occurred at the point (0 C ,O). We were 
able to determine the explicit form of the solution curve (16.31) by applying a Liapunov- 
Schmidt reduction to (16.21) at the point (0 C) 0), which was possible given that the operator 
D^G{4> C , 0) had a one-dimensional kernel. Similarly, in Section l5T2l we showed that when 
the mean curvature r was an arbitrary, continuously differentiable function, solutions to 
F(((j), w), A) = were non-unique. Again, this followed because we explicitly computed 
the solution curve through the point ((0 C , 0), 0). In Section |5T2l we found that the solution 
curve through ((0 C , 0), 0) had the form 

4>(s) = <p c + s + ^\(0)u(x)s 2 + 0(s 3 ), (6.4) 
w (s) = ^X(0)v(x)s 2 + O(s 3 ), (6.5) 

\( s ) = l\(0)s 2 + O(s 3 ), (6.6) 

which we demonstrated by applying a Liapunov-Schmidt reduction to the system (16.11) at 
the point ((0 C , 0), 0). Again, this was possible because of our work in Section |4~T1 where 
we showed that the linearization DxF(((p c , 0), 0) had a one-dimensional kernel. 

The importance of these non-uniqueness results is that they demonstrate first and fore- 
most that the conformal formulation with unsealed source terms is undesirable given 
that solutions for this formulation will not allow us to uniquely parametrize physical so- 
lutions to the Einstein constraint equations. Additionally, this paper helps build on the 
work of Walsh in [19 ] by expanding the understanding of how bifurcation techniques can 
be applied to the various conformal formulations of the constraint equations. This work 
is also interesting in that the analysis conducted here helps clarify the ideas of Baum- 
garte, O'Murchadha, and Pfeiffer in fl4) by showing how terms with "the wrong sign" 
that contribute to the non-monotonicity (non-convexity of the corresponding energy) of 
the nonlinearity in the Hamiltonian constraint directly contribute to the non-uniqueness 
of solutions. Finally, it is hope of the authors that this work will also help to lay the foun- 
dation for future analysis of the uniqueness properties of the Conformal Thin Sandwich 
method and the far-from-CMC solution framework established in [|9l[T0l. 

7. Appendix 

7.1. Banach Calculus and the Implicit Function Theorem. Here we give a brief re- 
view of some basic tools from functional analysis. The following results are presented 
without proof and are taken from [20] . We begin with some notation. 

Suppose that X and Y are Banach spaces and U C X is a neighborhood of 0. For a 
given map / : U C X — > Y, we say that 

f(x) = o(||x||), x -¥ iff r(x)/||x|| ->• as x — »• 0. 
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We write L(X, Y) for the class of continuous linear maps between the Banach spaces X 
mdY. 

Definition 7.1. Let U C X be a neighborhood of x and suppose that X and Y are 
Banach spaces. 

(1) We say that a map f : U — >■ Y is F-differentiable or Frechet differentiable at 

x iff there exists a map T G L(X, Y) such that 

f(x + h)- f(x) = Th + o(\\h\\), as h -> 0, 

for all h in some neighborhood of zero. If it exists, T is called the F-derivative or 
Frechet derivative of f and we define f'(x) = T. If f is Frechet differentiable 
for all x G U we say that f is Frechet differentiable in U. Finally, we define the 
F-differential at x to be df(x; h) = f'(x)h. 

(2) The map f is G-differentiable or Gateaux differentiable at x iff there exists a 
map T G L(X, Y) such that 

f(x + tk) - f(x) =tTk + o(t), ast^O, 

for all k with \\k\\ = 1 and all real numbers t in some neighborhood of zero. 
If it exists, T is called the G-derivative or Gateaux derivative of f and we 
define f'(x) = T. If f is G-differential for all x G U we say that f is Gateaux 
differentiable in U. The G-differential at x is defined to be dcf{x; h) = f'(x)h. 

Remark 7.2. Clearly if an operator is F-differentiable, then it must also be 
G-differentiable. Moreover, if the G-derivative f exists in some neighborhood of x and 
f is continuous at x, then f'(x) is also the F-derivative. This fact is quite useful for 
computing F -derivatives given that G-derivatives are easier to compute. See [|20l for a 
complete discussion. 

We view F-derivatives and G-derivatives as linear maps f'(x) : U — > L(X,Y). More 
generally, we may consider higher order derivatives maps of /. For example, the map 
f"(x) : U — > L(X, L(X, Y)) is a bilinear form. We now state some basic properties of 
F-derivatives. All of the following properties also hold for G-derivatives. 

The Frechet derivative satisfies many of the usual properties that we are accustomed 
to by doing calculus in IR n . For example, we have the chain rule. 

Proposition 7.3 (Chain Rule). Suppose that X, Y and Z are Banach spaces and assume 
that f : U C X — > Y and g : V C Y — )■ Z are differentiable on U and V resp. and 
that f(U) C V. Then the function H(x) — g o f, i.e. H(x) = g(f(x)), is differentiable 
where 

H'{x) = g'(f(x))f(x) 
where we write g'(f(x))f'(x) for g'(f(x)) o f'(x). 

Given an operator / : X x Y — > Z, we can also consider the partial derivative of / 
with respect to either x or y. If we fix the variable y and define g{x) = f(x, y) : X — >■ Z 
and g(x) is Frechet differentiable at x, then the partial derivative of / with respect to x 
at (x, y) is f x (x, y) = g'(x). We can a make a similar definition for f y (x, y). Finally, we 
observe that we can express the F-differential of f'(x, y) in terms of the partials by using 
the following formula: 

f'(x, y)(h, k) = f x (x, y)h + f y (x, y)k. (7.1) 

We have the following relationship between the partial derivatives and the Frechet 
derivative. 
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Proposition 7.4. Suppose that f : X x Y — >• Z is F-differentiable at (x, y). Then the 
partial F-derivatives f x and f y exist at (x, y) and they satisfy (17.11 ). Moreover, if f x and 
f y both exist and are continuous in a neighborhood of (x, y) then f'(x, y) exists as an 
F-derivative and (f7.ll) holds. 

7.1.1. Implicit Function Theorem. Suppose that F : U x V — > Z is a mapping with 
U C X, V C Y and X, Y, Z are real Banach spaces. The Implicit Function Theorem 
is an extremely important tool in analyzing the nonlinear problem 

F(x,y) = 0. (7.2) 

We present the statement of the Theorem here, the form of which is taken from [fT3l . For 
a proof see 11201161 

Theorem 7.5. Let (|7.2I) have a solution (xo, yo) £ U x V such that the Frechet derivative 
of F with respect to x at (x , y ) is bijective: 

F(x ,y ) = 0, (7.3) 

D x F(xq, yo) :— > Z is bounded (continuous) 

with bounded inverse. 

Assume also that F and D X F are continuous: 

F eC(U XV,Z), (7.4) 

D X F e C(U x V, L{X, Z)), where L(X, Z) 

denotes the Banach space of bounded linear operators 

from X into Z endowed with the operator norm. 

Then there is a neighborhood U\ x V\ C U x V of(x , y ) and a map f : V\ — > U% C X 
such that 

/(yo) = x , (7.5) 
F(f(y),y) = for ally eV v 

Furthermore, f E C(Vj . X) and every solution to (17.21) in U\ x V] is of the form ( f(y).y). 
Finally, if F is k-times differentiable, then f is k-times differentiable. 

7.2. Elliptic PDE tools. Here we assemble some useful tools for working with nonlin- 
ear elliptic partial differential equations. Throughout this section we will assume that Ai 
is a closed manifold with a smooth SPD metric g a b and that A is the associated Laplace- 
Beltrami operator. 

7.2. 1 . Maximum Principle. In this section we present a version of the maximum princi- 
ple on closed manifolds. The following result is well-known, but we present it here for 
completeness. 

Theorem 7.6. Let u e C 2 (M). Then if 

Au > or Au = or Au < 0, (7.6) 
then u must be a constant. In particular, the problem 

Au = f(x,u), 

has no solution iff(x, u) > or f(x, u) < unless f(x, u) = 0. 

Proof. See [fT8l for a proof. □ 
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7.2.2. Method of Sub- and Super-Solutions. Here we present a theorem that provides a 
method to solve an elliptic problem of the form 

Lu = f(x,u), (7.7) 

where 

Lu = -Au + c{x)u, c(x) e C{M x R) , c(x) > (7.8) 
and the function f(x, y) is nonlinear in the variable y. 

Theorem 7.7. Suppose that f : M x R+ ->f is in C k (M xR + ). Let L be of the form 
(17.81 ) and suppose that there exist functions w_ : M. — > R an J w + : .M — >■ R smc/? ?/j<3? 
the following hold: 

(1) u_,u + e C k {M), 

(2) < u-(x) < u+{x) Vx G A^, 
(5J Ln_ < f(x,u-), 

(4) Lu + > f(x,u+). 
Then there exists a solution u to 

Lu = f(x,u) onM, (7.9) 

such that 

(i) u e C\M), 

(ii) u-(x) < u(x) < u + (x). 

Proof. See ifTTTl for a proof. □ 



7.2.3. Fredholm Properties and Liapunov -Schmidt Decompositions for 
Elliptic Operators. In this appendix we discuss the Fredholm properties of linear ellip- 
tic operators on a closed manifold. We use these properties to form Liapunov-Schmidt 
decompositions for a given elliptic operator L between certain Banach spaces. The fol- 
lowing treatment is taken from [[TBI . 

Let u G C 2 ' a (M) and define the elliptic operator L : C 2 > a {M) ->■ C^ a {M) by 

n n 

Lu = - 22(aij{x)u Xi ) X j + / J bj{x)u Xi + c(x)u, (7.10) 

i,j=l i=l 

where a^ , 6j and c are smooth, bounded coefficients where = a^. We also assume 
that the satisfy the standard elliptic property 

n 

aiMi > d U\\\ 

where d > is constant and || • || is the Euclidean norm on R n . 
The operator (17.101 ) has an associated bilinear form 

B(u,u) = (Lu,u) = (u,L*u), (7.11) 

where (•, •) is the L 2 (M) inner product and L* is the L 2 -adjoint defined by 

n n 

L*u - -} (>/,;(./•)//,. - s ^2(b i {x)u) Xt + c{x)u. (7.12) 

i,j=X i=l 

Using the bilinear form B(u, u), the elliptic operator (17.101 ) defines an elliptic 
operator 

L : L 2 (M) L 2 (M), with domain of definition D(L) = H 2 (M). (7.13) 
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It is a standard argument in linear elliptic PDE to show that there exists a c > such the 
operator L + cl : H 2 (M) — >■ L 2 (Ai) is bounded and bijective. In particular, one shows 
that there exists a c > such that the associated bilinear form B(u, u)+c\\u\\2 is coercive 
and then applies the Lax-Milgram Theorem to conclude that there exists a unique weak 
solution u G H X {M) to 

Lu — cu = f for every / £ L 2 (Ai). 

Standard elliptic regularity theory implies that u & H 2 (M) and the norm || ■ || 2 ,2 makes 
D(L) a Hilbert space. An application of the Open Mapping Theorem (Bounded Inverse 
Theorem) then implies that 

(L + ciy 1 : L 2 (M) -> D(L), 

is continuous. This implies that the operator L + cl is closed and that the operator 
(L + cl) — cl = L is closed. In addition, the operator 

K C = (L + cl)' 1 G L(L 2 (M), L 2 (M)) is compact 

given that the embedding H 2 (M) C L 2 (.M) is compact. For / £ L 2 (M), we have the 
equivalence 

Lu = f, ueH 2 (M)^ (7.14) 
u - ciT c M = KJ, u e L 2 (M). (7.15) 

Riesz-Schauder theory implies that (/ — cK c ) is a Fredholm operator and the equivalence 
(17.141 ) implies that L is a Fredholm operator. 

Because L is a Fredholm operator of index zero, we have that R(L) is closed. There- 
fore we may write 

L 2 (M) = R{L) © Z Q , 

where Z = R(L) 1 - is the orthogonal complement with respect to the LMnner product. 
Because D(L) is dense in L 2 (M) and L is closed, may apply the Closed Range Theorem 
to conclude that 

R(L) = {f e L 2 (M) | (/, u) = for all u e N(L*)} (7.16) 
and that Z = N(L*), where L* : L 2 (M) -)> L 2 (>1) is induced by (I7T21) . Therefore 

L 2 (7W) = J R(L)©iV(L*), 

and if D(L*) = H 2 (Ai), the above arguments imply that L* is Fredholm operator. So 
we have the following decomposition of the codomain of L*\ 

L 2 (M) = R(L*)®N(L). (7.17) 

Finally, given that N(L) C D(L) = H 2 (M) C L 2 (M), the decomposition (17.171) al- 
lows us to obtain the following Liapunov-Schmidt decomposition for the linear problem 

L:H 2 (M)^L 2 (M): 

H 2 {M) = N(L)®(R(L*)nH 2 (M)), (7.18) 
L 2 (M) = R{L) ®N{L*). (7.19) 

Now we observe that the Fredholm properties of linear elliptic operators derived on 
Hilbert spaces hold for subspaces that are only Banach spaces. We then use these Fred- 
holm properties to derive Liapunov-Schmidt decompositions for these Banach spaces. 
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Suppose that the Banach space Z C L 2 (Ai) is continuously embedded and that the 
domain of definition X C Z with a given norm is a Banach space that satisfies the 
following conditions: 

L : X — > Z is continuous , (7.20) 

Lu = f foru E D(L) = H 2 (M)J E Z E X. 

Equation (|7.20l) is an elliptic regularity condition and is satisfied for a variety of spaces, 
most notably X = W 2 ' P (M), Z = U\M) and X = C 2 ' a (M), Z = C°' a (M) with the 
standard norms. Then for X and Z satisfying (17.181 ) and (17.201 ) we have that 

N(L) = N(L\z) C X, and (7.21) 
R(L) n Z = R(L\ Z ) is closed in Z, (7.22) 

given that Z C L 2 (Ai) is continuously embedded and R(L) is closed in L 2 (M). The 
ellipticity property (17.201) also holds for the adjoint L* and implies that 

N(L*) C X, where D(L*) = D(L) = X. 

Applying the decomposition (17.191) , we may write any z E Z as 

z = Lu + u*, where u G D(L),u* E N(L*), (7.23) 

Lu = z — u* E Z =^> u E X, therefore 

Z = R{L\ Z )®N{L*). 

Finally, we have that dimiV(L| z ) = dimiV(L) = dimiV(L*) and that 

L : X — > Z, X = D(L\z), is a Fredholm operator of index zero. (7.24) 

The decomposition (17.181) then implies that 

X = N{L\ Z ) © (R(L*) n X), (7.25) 

and so (17.231) and (17.251) constitute a Liapunov-Schmidt decomposition of the spaces X 
and Z with respect to a given linear, elliptic operator L. 

Remark 7.8. As noted in H3]|, we may regard the spaces W 2 ' P (M) C L P (M) C L 2 (M) 
for p > 2, and we can then apply the above discussion to conclude that a linear elliptic 
operator L : W 2)P (M) — > L p (Ai) is Fredholm and use this fact to obtain a Liapunov- 
Schmidt decomposition of X = W 2 ' P (M) and Z = L P (M). Similarly, C 2 > a (M) C 
C°> a (M) C L 2 (M)fora E (0,1), so L : C 2 > a (M) -)■ C°> a (M) is Fredholm and 
we may also obtain a Liapunov-Schmidt decomposition of X = C 2,a (M) and Z = 
C°> a (M ) using (17331 and (17351) . 
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